Bifurcation Phenomena in Optimal Velocity Model for Traffic Flows 
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In the optimal velocity model with a time lag, we show that there appear multiple exact solutions 
in some ranges of car density, describing a uniform flow, a stable and an unstable congested flows. 
This establishes the presence of subcritical Hopf bifurcations. Our analytic results have far-reaching 
implications for traffic flows such as hysteresis phenomena associated with discontinuous transitions 
between uniform and congested flows. 
PACS numbers: 05.45.Yv, 04.20.Jb, 45.70.Vn, 47.54.+r 



The optimal velocity (OV) model p| for traffic flows 
attracts considerable interest recently, since it describes 
congested flows as well as uniform flows. For low mean 
densities, the uniform flow is realized, while it becomes 
unstable above some critical density and the congested 
flow shows up, with high and low density regions ap- 
pearing alternately. At the critical density, the model 
predicts the presence of a discontinuity in the correla- 
tion diagram between car density and traffic flux, called 
the fundamental diagram. This discontinuous behavior 
is consistent with observed data. (In Fig. 9 of their 
results are compared to an observed data ||.) 

Although it is an important consequence of the OV 
model, the precise origin of the discontinuity has not been 
identified. The purpose of this letter is to demonstrate 
that the subcritical Hopf bifurcation phenomena Q is 
the dynamical origin of the discontinuous transitions be- 
tween uniform and congested flows. We consider in this 
letter an OV model with a time lag whose exact solutions 
have been found recently by three of the present authors 
Q(see also ref. Q). We observe that in some ranges of 
density three exact solutions coexist; a stable uniform 
flow, a stable and an unstable congested flows. In a nu- 
merical simulation one of them is realized depending on 
its initial condition. Any congested flow, where all cars 
have the same maximum (u ma x) as well as the minimum 
value of velocities (v m i a ), may be characterized by the 
amplitude in velocity, Aw = v max — u m in- As for the coex- 
isting solutions, Av of the stable congested flow is larger 
than that of the unstable flow, while Av — for the uni- 
form flow. These values are functions of mean headway 
h and can be used to establish the bifurcation to be of 
subcritical type. Further the discontinuity in the traf- 
fic flux mentioned earlier is interpreted as the hysteresis 



phenomenon associated with the subcritical bifurcation. 
The coexistence of stable uniform and congested flows 
H and the hysteresis phenomena |7]] were noticed ear- 
lier in simulations. However let us emphasize here that 
the subcritical bifurcation is the key concept to explain 
properties of the transition between the two flows. 

The system of first-order differential-difference equa- 
tions H Q we consider is given by 



i n {t + t) = V[Ax n (t)] = £ + nt&nh 
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where x n and Ax n = x n ^\ — x n correspond to the posi- 
tion of nth car and its headway, the distance between the 
car and the preceding (n — l)th car. r is the time lag to 
reach the optimal velocity V(Ax) when the traffic flow 
changes. We impose the periodic boundary conditions, 
x n +N = x n — L, where N is the total number of cars in a 
circuit of length L. The OV function V(Ax) is described 
by a hyperbolic tangent with four parameters, £, rj, p and 
a. 

The trivial solution of the system (||) corresponding to 
a uniform flow is given by 
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(t) = V(h) t-nh, 
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where h = L/N is the common headway. A linear stabil- 
ity analysis shows that the uniform flow is unstable for 
the region | h — p \ < Hq. The critical values p ± H are 
determined by M 
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where r c = a/rj. Thus, for — H + p < h < H + p, an 
initially prepared uniform flow becomes unstable, and a 
congested flow develops. As shown in ref. ||, the result- 
ing congested flow is described by an analytic function of 
the form 



x n (t) =Ct -nh + A In 
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where A, (3, v, C, 5 are constants. $q(v, q) is one of the 
theta functions with the modulus parameter q. It follows 
from the periodicity that 2(5 = rib/N, where an integer 
rii is the number of low density regions in the circuit. 
The width parameter < 2<5 < 1 determines the ratio 
of the length of a low density region to L/nf. For each 
length L/rib, the length of the low density region is given 
by 25 L/ rib, while that of the high density region by (1 — 
28)L/rib. In this sense, by replacing 26 by 1 — 25, we 
exchange the length of two regions. 

A crucial point for observing bifurcation is the recogni- 
tion of rich structure of the exact solution (^). The stable 
uniform flows always appear in the region H < | h — p \ . 
Unlike this, there must be critical values of headway 
p ± Hi (Hi > H ) such that the periodic solutions do 
not exist for h — p |> Hi: no congested flows are 
generated if the mean density is too high or too low. 
Thus, h can be divided into five regions: h > Hi + p (I), 
H a + p<h<Hi+p (II), -H + p<h<H a + p (III), 
-Hi + p < h < -H + p (IV) and -Hi +p> h (V). 
We have uniform flow solutions over the whole region; 
They are stable except in the region (III). There are no 
congested flow solutions in (I) and (V); The stable ones 
are known to appear in (III). Since the stability of the 
uniform flows in (II) and (IV) is confirmed against small 
perturbations only, one cannot exclude the presence of 
other stable solutions than the uniform flows. We shall 
show that stable as well as unstable periodic solutions 
do coexist actually in these regions, and that both are 
described by (|]). 

To this end, we discuss the relations of the model pa- 
rameters (r, £, rj, p, <j, N, L) in (Q) to those in the ansatz 
(§), (A,vJ3,8,q,C,h). In addition to L = Nh the fol- 
lowing |l7J] should obey in order for (^) to solve (|l|); 
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Inclusion of the relation 2/3 = rib/N with a given rib 
makes seven reltaions between the above two sets of seven 
parameters |10 : thus, we may construct exact solutions 
for a given set of the model parameters. In this paper 



we restrict ourselves to the case rib — 1, f° r simplicity. 
(See ref. [jlfj for solutions with lib / 1.) The eq.@ is the 
Whitham's dispersion relation |J, an important charac- 
teristic of the delayed model. Note that h — p and C — £ 
change their signs while rj remains invariant when we re- 
place 25 by 1 - 28. 

We show the presence of multiple solutions for a given 
set, (r, £, i], p, a, N, L). To obtain the parameters 
of the solutions, one first fixes (A, v, P) from (||), (||) 
and 2/3 N — 1. As the modulus parameter, we use k = 
— 7r/(ln<7) instead of q. By solving (^), 5 is computed as a 
function of n. (As a result of the property of 5 discussed 
earlier, it has two branches, 5i(k) and 62 (k) such that 
25i(k) = 1 — 2$2(k).) Then substituting 5(k) to (||) gives 
us a function, /i(re). The result is shown in Fig.l. 
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FIG. 1. The relation between the mean headway h and 
ft. The solid and broken lines correspond to stable and 
unstable solutions, respectively. Here the parameters are 
t = 0.58228, £ = tanh2, 77 = 1, p = 2, a = 1/2 and N = 20, 
for which H = 0.38978 and Hi = 0.56290. 

As shown in Fig.l, its graph is symmetric around the 
axis h — p. At the symmetric point, 25 becomes 1/2. It 
should be noted that n(h) becomes a two- valued function 
in the regions (II) and (IV). In each region, we have two 
periodic solutions for a given set of the model parame- 
ters. We have confirmed by numerical simulations that 
one (denoted by A on the solid line in Fig.l) is stable, 
while the other (denoted by B on the broken line in Fig.l) 
is unstable: When a small perturbation is added, the lat- 
ter goes down to the former or to a uniform flow with the 
same value of h. Therefore, there are altogether three so- 
lutions in this region, indicating a subcritical bifurcation. 
The bifurcation points h = p ± Ho arc determined from 
(||) and ^ for k = 0, which end up with (||). The value of 
k which gives the critical values h = p±Hi is determined 
by 
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where B = 2_ftT/3, Z is the Jacobi's zeta function and K 
is the complete elliptic integral of the first kind. This 
determines Hi. 

As an "order parameter" characterizing the subcritical 
bifurcations, we may take the velocity amplitudes Av of 
relevant flows as functions of h. In Fig. 2 we have drawn 
the analytical results with the solid and broken lines, 
which agree with numerical simulations shown with the 
small squares. The simulations are done by gradually 
increasing h and the results for Av change as indicated 
with thin arrows: when we increase h, there is a jump in 
Av at h = p—Ho, from a uniform flow to a congested flow; 
at h = p+Hi there is another jump from a congested flow 
to a uniform flow. (Although not indicated in the figure, 
jumps occur at h = p + Hq and h = p — Hi for decreasing 
h.) This is the hysteresis phenomenon associated with 
the subcritical bifurcation. It is appropriate to explain 
how we performed the simulations. Suppose we have a 
stable flow as a result of the simulation at a value for 
the mean headway h, say ho- Taking a snapshot of the 
flow, we know the positions of cars at the time. Then we 
measure the headway of each car and increase it by Ah, 
and that configuration is used as the initial condition for 
the simulation with h = ho + Ah. This procedure gives 
the results for the region (II) staying on the curve for 
the congested flow: it is important to keep the Ah small 
enough since the results would heavily depend on the 
initially prepared configurations |j. 
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FIG. 2. The h-Av relation from the exact solution, com- 
pared with a numerical study. The simulations were done by 
increasing h as indicated by arrows, the results are plotted 
with the small squares. We observe a hysteresis phenomenon. 



We discuss now flux-density correlation diagram, the 
fundamental diagram. The flux Q may be defined by the 
number of cars passing through some reference point in 
a certain time interval. For a uniform flow, it is given 
by product of the density 1/h and the common velocity 
V{h): 



Q = V(h)/h. 



(11) 



We calculate the flux of a congested flow as follows. Since 
high density regions are moving backward with velocity 
vb, it will be convenient to take the rest frame of the 
density waves. The period To defined for a car to make 
a round in the rest frame is given by 



Tn = - = 2r N, 



(12) 



where we have used the Whitham's relation, (3 = tv — 
1/(2 A). Since the average velocity of cars in the original 
coordinate is to be L/Tq — vb, the flux Q may be given 
by 



\{y - VB ) = h~ VB,h - 



(13) 



Let us calculate vb by considering the repetitive spatio- 
temporal pattern in a congested flow: 



i(t) = x n (t + 2r) + 2 v B t. 



(14) 



This implies that the x n has the same time dependence 
as x n —i apart from a time delay 2r and a distance 2vbt. 
It follows from (g) and (pi) that 



Therefore, one obtains a simple expression for Q 
Q = C/h. 



(15) 



(16) 



As discussed previously, one calculates the parameter C 
as a function of k using (R). Together with h(k), it gives 
C(h). Our analytic results are shown in Fig. 3: Disconti- 
nuities in the traffic flow show up as a result of a hystere- 
sis phenomenon, which is consistent with the observed 
data qualitatively. Since the size of the discontinuous 
jumps may be determined analytically, a quantitive com- 
parison is also possible. 

A few comments are in order. (1) The presence of the 
subcritical Hopf bifurcations naturally define the bound- 
aries between (I) and (II), and between (IV) and (V). 
This should be compared with the Fig. 8 of (2) In 
Fig. 3, the discontinuities in the low density side appears 
after passing through the uniform-flow's peak. Since the 
position of the peak is determined by the OV function, 
we may take a set of model parameters so that the dis- 
continuities appear before the peak. 
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FIG. 3. The fundamental diagram for traffic flow. The 
dotted line is for uniform flows, while the solid and broken 
lines are for stable and unstable congested flows, respectively. 
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In summary, the presence of subcritical Hopf bifurca- 
tions in a delayed optimal velocity model which admits 
exact solutions is established by showing coexistence of 
stable and unstable solutions. Although our results are 
obtained in a specific OV model, we believe that any OV 
type model will show up the subcritical bifurcations. Ac- 
tually, we have confirmed by numerical simulations the 
presence of multiple stable solutions for the model given 
by [V(Ax n ) - in] (First studied in in this con- 

text). Therefore, our results given here should be con- 
sidered, at least qualitatively, to be a universal feature 
of the OV models describing spatio-temporal patterns. 
They may be of fundamental importance in the area of 
the traffic engineering. In order to confirm the subcriti- 
cal bifurcation directly highly desirable are more precise 
data of the h-Av relation. 

Parts of our numerical works were done with the 
computer fascilities of Yukawa Institute for Theoretical 
Physics. We are grateful to K. Nishinari, H. Hayakawa 
and Y. Ono for discussions and suggestions. 
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